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The closed form of the first order non-linear differential equation that is satisfied by the effective 
range within the variable phase formulation of scattering theory is discussed, ft is shown that 
the conventional method of determining the effective range, by fitting a numerical solution of the 
Schrodinger equation to known asymptotic boundary conditions, can be modified to include the first 
order contribution of a long range interaction. 



I. INTRODUCTION 



Scattering of atoms and molecules in an ultra-cold environment is a process that is very important in the interpreta- 
tion of measurements of properties of cold trapped gases and in the study of Bose-Einstein condensation. At such low 
temperatures s-wave scattering dominates and the scattering length, a, usually suffices to describe collisions. However 
the effective range, r e , is important because it implies that there is a relative contribution to the elastic scattering 
cross section at energy E of 2fiEar c where fx is the reduced mass. The effective range can provide an upper bound to 
the temperature, T, of a trapped sample in thermal equilibrium below which only the scattering length is significant; 
for example the corresponding relative contribution to the elastic scattering cross section that is attributable to the 
effective range is 2iik^Tar c where fee is Boltzmann's constant. Previously [l| we used the variable phase formulation 
of scattering theory [2 to study the calculation of the scattering length in the presence of long range forces. Here we 
make a study of the effective range. 

In the conventional method of calculating the effective range a numerical solution of the Schrodinger equation for 
scattering is fitted to known asymptotic boundary conditions at some suitably large separation i? c to yield a phase 
shift. This is repeated for several small values of the asymptotic wave number of the relative motion, k, and the 
values of a and r e are obtained from the effective range theory fit of the phase shift as a function of k 0. This 
procedure is appropriate even in the presence of long range forces because, by making the fit at the finite separation 
R c , we solve a problem for short ranges forces curtailed at R c - A difficulty occurs in that R c must be large. Here 
we present a correction that accounts for the influence of a long range force over the range [R Cl oo[ to first order in 
the potential strength. It allows much smaller values of R c to be used with attendant reduction in computation time 
and accumulated error. Marinescu y| has provided first and higher order corrections but they are complicated in 
comparison to our formulation and, in most cases, a first order correction is enough. 

II. THEORY 

In the notation of Levy and Keller |5j the s-wave phase shift 77 is expanded as 

ri = Ntt + k(ri + fc 2 r/i + . . .) (1) 

where N is the number of bound states supported by the interaction potential. Equation (TlJ is in accord with 
Levinson's theorem. Levy and Keller [f| omitted the term Ntt from their equation (8) but the subsequent analysis is 
not altered because effective range theory 3] provides an expansion for k cot r\. 

When long range forces are present expansion and the effective range expansion have additional terms in k, 
such as fc 4 In k, but for interactions dominated by a Van der Waals R~ 6 dispersion potential the leading terms have 
the indicated fc-dependence although the coefficients, T)i, are modified jfj]. In variable phase theory the phase shift 
and these coefficients are functions of the separation, R, and hence the scattering length and the effective range are 
also functions of R. The effective range expansion is 

kcot V (R) = --^ + ^r e (R)k 2 + .... (2) 

The scattering length, a(R), is the limit as k — > of — ta,nrj(R) / k and hence is —r)o(R). The effective range is twice 
the coefficient of the term in k 2 in equation J5J and from the Taylor expansion of the cotangent it can be expressed as 

2a(R) 2 Vl (R) 

r c (R) = — - wfip- (3) 
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The usual scattering length and effective range are the asymptotic values taken by a(R) and r (R) as R — * oo. 
The log-derivative in variable phase theory is 



u{R) 



dr)(R) 



dR 



cot[kR + r)(R)]. 



(4) 



Substitution of u(R) into the Riccati equation that is satisfied by the log-derivative yields the variable phase 
equation (equation (6) of Levy and Keller Q) 



MR) 

dR 



-V(R) sm 2 [kR + i](R)} 



(5) 



where V(R) denotes 2/i/h 2 times the interaction potential. From equation 0) Levy and Keller [§J rearranged equation 
(J5J as a power series in k to obtain non-linear first order differential equations in R that are satisfied by the coefficients 
r)i(R) (their equations (9) and (10)) from which we find 



da(R) 
dR 



and 



= V(R)[R-a{R)} 2 



(6) 



= -2V(R)[R - a(R)MR) + ^V(R)[R - <ffi 4 - 



(7) 



A. Numerical procedures 

Although equation Q has the closed form solution, 



m(R) 



exp 



V(S)[S-a(S)]dS } x 



exp ■ 



V(Q) [Q - a(Q)}dQ } V(S) [S - a(S)] 4 dS, (8) 



we find that numerical evaluation is not conveniently provided in this way because the closed form solution (jSJ requires 
computation of a double integral and, when the interaction supports bound states, also requires a suitable numerical 
account of the poles in the scattering length, a(R); these poles also make a direct numerical solution of equation @ 
impractical. The adaptive method of Lambert Q is unsuitable for heavy particle collisions because there are many 
poles in a(R). The log-derivative method that we used previously [l| is also unsuitable; although equations © and 
may be written in the Riccati form 



d 

dR 



R-a(R) r]i(R) 
R — a(R) 



V(R) 



R — a(R) ?7i(i?) 
R - a(R) 



1 ^V{R)[R-a(R)} 4 
1 



(9) 



the difference between the right hand side and the unit matrix prevents application of the log-derivative method. A 
rational fraction transformation of equation @ yields another Riccati equation [8j but there is no obvious transfor- 
mation that simultaneously changes the right hand side of equation @ to a unit matrix and maintains the structure 
of its left hand side. 

We suggest that the conventional method for obtaining the effective range is used at separation R Cl chosen beyond 
the last pole in a(R), and equations and Q) are then used to obtain a correction that is of first order in the 
potential strength. 



III. LONG RANGE CORRECTIONS 



The exponential terms in expression (jSJ are positive. Hence the contribution to rji(R) from the range [R Cl oo) is 
negative when the potential is attractive over this range and therefore the conventional calculation yields an upper 
bound to ?7i = 771(00). 
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The subsequent analysis is simpler when presented in terms of dimensionless quantities obtained by scaling by 
suitable powers of R c . Throughout the remainder of this note, unless otherwise indicated, we let R, a, r c , V and t]i 
denote the corresponding quantities divided by R c , R c , R c , R~ 2 and R^ respectively. Equations J3J), JEJ and JJJ) are 
not altered. The correction, Sr]i, to rji, obtained by integrating the right hand side of equation (JJ by parts over the 
range [1, 00 [, is 



5th = -%Vi \v w (a - 1) + V {2 A 



4! 



where 



V {i) (R) 



-(5-i) 



V {l ~ 1) {S)dS 



■V (l \a-1) 



5-t 



(10) 



(11) 



with i = 1,2, ... ,5 and V^(R) identified as V(R), and the right hand expression is evaluated at R — 1 (i.e. at 
separation R c ). 

The inverse power potential that is given in terms of the unsealed separation by —C n R~ n is represented in the 
scaled quantities by 



V(R) = % 



where 



( 2nC n 



R" 



,-2 1 
Rc. 



(12) 



(13) 



From equations ©, ©■ <|TU|) . (fTTft . (|T2*)l and we obtain the first order correction, 5r e , to the effective range as 



5r e 



„n— 2 



2r„ 



2r fi 



(n — 3)a (n — 2) 3(n — 5)a 2 3(n — 4)a n — 3 



(14) 



where the right hand side is evaluated at R = 1. Equation i|14|) is consistent with the correction that was derived by 
different methods [|j • Note that the correction given by equation (|14H is not identical to the correction in reference ^ ; 
equation l|14|) provides the correction to the effective range whereas reference gives the correction to the product 
r c a 2 . Their equivalence is readily shown by introducing the correction to the scattering length; the correction is 9] 



Sa 



,71-2 



1 



2 a 



^2 1 



n — 3 n — 2 n — 1 



(15) 



IV. METHOD 



Our suggested procedure for calculating the effective range is as follows. First evaluate the scattering length and 
effective range at separation R c by the conventional method and then scale them by division by R c . For each term 
such as (|12H of the long range potential, evaluate the scaled correction 5r c of equation i|14|) and find the sum of 
the corrections; this is legitimate because we are interested in only the first order correction. Multiply the resulting 
correction by R c . Some experiment is needed to find a suitable value for R c but generally it is much smaller than the 
separation that represents infinity in the conventional method without correction. The leading terms that describe 
the influence of the long range force that are not accounted for by equation H14[) are of second order and hence are 
proportional to (a™~ 2 ) 2 . The chosen value of R c should coerce this quantity to be smaller than unity; a suitable 

choice is several times the characteristic length, (2/zCn/fi 2 ) 1 ^" 2 , of the potential. 
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TABLE I: Effective Range and Corrections (bohr) 



Separation R c 
(bohr) 


Effective Range 


c 6 


Correction from term 
C 8 


in 

Cio 


Effective Range 
(Corrected) 


600 


284 


317 


428 xlO" 4 


1088 xl0~ B 


602 


800 


358 


259 


200 xl0~ 4 


287 xlO" 8 


617 


1000 


406 


215 


107 xl0~ 4 


99 xl0~ 8 


621 


1200 


440 


183 


64 xlO" 4 


41 xl0~ 8 


623 


1400 


464 


160 


41 xlO" 4 


19 xl0~ 8 


624 


1600 


483 


141 


28 xlO" 4 


10 xlO" 8 


624 


1800 


498 


126 


20 xlO" 4 


6 xl0~ 8 


624 


2000 


510 


115 


15 xlO" 4 


3 xlO" 8 


624 


3000 


546 


78 


4 xl0~ 4 


< 10~ 8 


624 


4000 


565 


59 


2 xl0~ 4 


< 10" 8 


624 


5000 


576 


47 


10" 4 


< 10~ 8 


624 


6000 


584 


40 


< 10" 4 


< 10 -8 


624 



V. EXAMPLE 

For the model Cs-Cs potential described by Gribakin and Flambaum |l(J we obtained the effective range accurate 
to three significant figures with R c = 1400 bohr while even the first significant figure is uncertain in the conventional 
uncorrected calculation with R c = 10000 bohr; the characteristic length of the dominant term of the asymptotic 
dispersion potential is 203 bohr. We show details of the convergence in tabled the values shown in tableware unsealed. 
The characteristic lengths of the dispersion terms in Cq, Cg and Cio are, in bohr, 203, 80 and 50 respectively. 

VI. CONCLUSION 

The variable phase method provides a simple way in which account may be taken of the long range dispersion terms 
of an interatomic potential in calculations of the effective range. The method gives considerable saving of computation 
time and of accumulated truncation error. 
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